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ABSTRACT 

Isolated compact groups of galaxies (CGs) present a range of dynamical states, group 
velocity dispersions, and galaxy morphologies with which to study galaxy evolution, 
particularly the properties of gas both within the galaxies and in the intragroup medium. 
As part of a large, multiwavelength examination of CGs, we present an archival study 
of diffuse X-ray emission in a subset of nine Hickson compact groups observed with the 
Chandra X-ray Observatory. We find that seven of the groups in our sample exhibit 
detectable diffuse emission. However, unlike large-scale emission in galaxy clusters, 
the diffuse features in the majority of the detected groups are linked to the individual 
galaxies, in the form of both plumes and halos likely as a result of star formation or 
AGN activity, as well as in emission from tidal features. Unlike previous studies from 
earlier X-ray missions, HCGs 31, 42, 59, and 92 are found to be consistent with the 
Lx~T relationship from clusters within the errors, while HCGs 16 and 31 are consistent 
with the cluster Lx — a relation, though this is likely coincidental given that the hot gas 
in these two systems is largely due to star formation. We find that Lx increases with 
decreasing group H I to dynamical-mass ratio with tentative evidence for a dependance 
in X-ray luminosity on H I morphology whereby systems with intragroup H I indicative 
of strong interactions are considerably more X-ray luminous than passively evolving 
groups. We also find a gap in the Lx of groups as a function of the total group specific 
star formation rate. Our findings suggest that the hot gas in these groups is not in 
hydrostatic equilibrium and these systems are not low-mass analogs of rich groups or 
clusters, with the possible exception of HCG 62. 
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1. Introduction 



As the majority of galaxies in the nearby Universe are found within gravitationally bound 
groups (e.g., Tully 1987; Small et al. 1999; Kar achentsev|2005 ) , understanding the physical processes 
at work in these systems is fundamental to galaxy formation, evolution, and cosmological theory. 
Though a significant fraction may be condensed structures embedded within largely extended, 



loosely bound systems (Tovmassian et al. 2006, Mendel et al. 2011), Hickson compact groups 



(HCGs; Hickson||l982 ) in particular are useful tools for studying tidally enhanced galaxy evolution 
in the nearby Universe because they are isolated from other nearby galaxies due to their selection 
criteria, have low velocity dispersions (usually < 200 km s _1 ), and exhibit high number densities 
(typically 3-4 galaxies of comparable luminosity within several galaxy radii of one another). These 
properties combined encourage numerous gravitational interactions between group members. 



Because the crossing time of CGs is short (< 0.02i# ; Diaferio et al.|1994 ), simulations indicate 
that group members should have entirely merged into single massive galaxies on timescales of 
~1 Gyr. Hypotheses concerning the continued existence of CGs have included that these systems 
are either recently formed (e.g., Mendes de Oliveira et al. 2005), that they are chance alignments 



within larger structures (e.g., Mamon||l986 ), or that CGs are perpetually formed within collapsing 



rich groups (e.g., Diaferio et al. 1994). However, all three proposed scenarios have evidence that 



contradicts them such as the wide range in ages of star clusters in CGs (Konstantopoulos et al. 



2010, 2012), the low surface brightness tidal debris and disturbed disk morphologies present in 
many groups (e.g., HCGs 16, 31, and 92; iMendes de Oliveira et al. 1119981 IGallagher et al. 112010 



Fedotov et al.|[20TT ), and indications that a small fraction of CG members are the result of mergers 
(Mendes de Oliveira et al. 2005). Clearly these systems are much more complicated than first 



suspected, and understanding the origin and physical processes responsible for maintaining CGs 
over timescales in excess of several Gyr is important for refining cosmological simulations. 



Verdes- Montenegro et al. (2001) presented an evolutionary scenario in which the spiral-rich 



compact groups with most of their neutral H I gas confined to the galaxy disks experience successive 
tidal encounters which liberate the cool gas from the galaxies into the intragroup medium in the 
form of tails and bridges. This progresses to one of two evolutionary end points: (1) a group in 
which the liberated H I is shock-heated to X-ray temperatures; or (2) a group with a large, low 
velocity dispersion H I halo. In a study of five groups, Freeland et al. ((2009) found that the amount 
of H I tidally removed from galaxy disks is related to the compactness of the group. Therefore, 
compact groups should be very efficient at dispersing their neutral gas into the intragroup medium. 
If a sufficient mass of gas is removed from the disks and heated to sufficiently temperatures, the 
group may be considered a fossil group. Jones et al. (2003) defined fossil groups as those with 



Lx > 10 42 h 5 Q erg s 1 and a difference in .R-band magnitude of Am > 2.0 mag between the two 
brightest group members. The possible evolution of HCGs into fossil groups represents the most 
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similar to clusters these systems may become. We note that tidal stripping of gas from galaxy 
disks is not necessarily the only source of baryons for the intragroup medium. Simulations of 
galaxy clusters at high redshift have shown that cold mode accretion of gas from the surrounding 
intergalactic medium can be a substantial source of baryons (e.g., Keres et al.||2005 ); however, it is 
unclear how significant this accretion may be for low redshift groups. 

Observations of diffuse X-ray emission are particularly helpful in placing constraints on the 
properties of CGs because the increased likelihood of tidal encounters implies that a significant 
fraction of the gas mass in these systems may be in the form of a hot plasma. Previous ROSAT 



(e-g., 


Ponman et al. 


1996. 


Mulchaey k Zabludoff 


1998 


Helsdon et al. 


2001, 


Mulchaey et al. 


2003) 


and Chandra (e.g., Rasmussen et al. 


2008, 


Rasmussen & Ponman 


2009 Sun et al. 


2009 


) analyses 



of diffuse X-ray emission in galaxy groups have specifically treated groups of galaxies as low- 
mass analogs of galaxy clusters. Early studies found that, within large errors, galaxy groups were 
consistent with the established relationships between bolometric X-ray luminosity Lx and gas 
temperature (Lx — T) as well as Lx and velocity dispersion (Lx — cr) found from observations 
of large samples of galaxy clusters. With improved instrumentation and more statistically robust 
datasets, it was shown that compact groups deviate from the established Lx — T cluster relation 
(groups are fainter for a given temperature), however the Lx — cr relation was still found to be 
consistent, though with a large scatter (cf. Figure 12 in |Ponman et al.||1996fr . 



Having performed the first extensive X-ray survey of compact groups using data from both 



the R OSA T All-Sky Survey (RASS) and pointed ROSATPSPC observations, |Ponman et al] (1996 > 
found that 22 out of 85 observed HCGs had diffuse X-ray emission above their detection limits. From 
the data, the authors statistically estimated that ~75% of HCGs have a diffuse X-ray luminosity 
above 1.3 x 10 41 erg s^ 1 . While diffuse X-ray emission was previously thought to be limited 



to groups containing only E/S0 members, Ponman et al. (1996) found that groups containing 



spiral galaxies exhibit diffuse X-ray emission as well, however the X-ray luminosity functions of 
spiral dominated groups show that they tend to be fainter than E/S0 dominated groups to a high 
statistical significance. Additionally, a weak anti-correlation was found between Lx and spiral 
galaxy fraction. 

Studies of compact groups such as HCG 31, which shows multiple interactions among its 
low-mass members, have shown evidence for hierarchical structure formation characteristic of en- 



vironments at higher redshifts (z~l-2) (Gallagher et al. 2010). Additionally, a gap in both the 



mid-IR colors and specific star formation rates (SSFRs; i.e., SFR normalized by stellar mass) of 
HCGs compared to galaxies in other environments suggests that the galaxies in these systems un- 
dergo accelerated evolution from star-forming "blue-cloud" galaxies to the quiescent "red sequence" 
d Johnson et aL|2007[ |Walker et aL]|2010[ |2012| |Tzanavaris et al.||2010[ ) . This rapid evolution is sup- 
ported by the H I deficiency observed in most HCG members relative to spiral galaxies in the field 



(Huchtmeier 1997 Verdes-Montenegro et al. 2001). 



The rapid evolution of HCG galaxies from gas-rich and star-forming to gas-poor and quiescent 
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raises the following question: how is the neutral gas being processed in HCGs, i.e., is most of the 
H I mass converted into stellar mass or is it ionized and expelled into the intragroup medium, 
and in what proportions do these mechanisms operate in individual groups? Previous papers have 
examined the star cluster populations in several HCGs (e.g., Gallagher et al.|2001 Konstantopoulos 



et al.pOlO Fedotov et al.poTT ) to study the consumption of gas by star formation, however there 



has been no investigation of how the H I gas is heated to X-ray temperatures in these systems 
in the context of group evolution. In this work, as part of a multiwavelength campaign to study 
the evolution of both the galaxies in compact groups and the group environment itself, we present 
the results of a study of the diffuse X-ray emission from nearby CGs (z < 0.015) using Chandra 
observations. Specifically, we inspect the morphology of the hot gas in compact groups using 
improved spatial resolution compared to previous studies that analyzed R OS A T data; we re-examine 
the Lx scaling relations with temperature and velocity dispersion; and we compare the hot gas in 
groups to the neutral H I gas to investigate how the X-ray gas is built up in these systems. 

First, we describe the Chandra observations and methods for reducing the events files in ^2} In 
we discuss the extraction of the X-ray spectra, as well as the models we used to fit the data. Our 
model-derived measurements are presented in §|4] in addition to a discussion of the results of our 
study in the larger context of galaxy evolution and the evolution of the group environment itself. 
We summarize our findings and discuss the future steps of our investigation into the diffuse X-ray 
properties of HCGs in Appendix [A] gives specific information on the extraction and modeling 
of the diffuse emission in the detected groups. Throughout this work, we assume a cosmology of 
Hq = 70 km s _1 Mpc -1 , Qyi = 0.27, and = 0.73. We also caution the reader that we allude 
several times to compact groups possibly becoming more cluster-like as they evolve. In this context, 
we refer to the hot gas distribution, not necessarily to the richness of the group; therefore, the term 
cluster-like is used to describe groups with a large fraction of E/S0 galaxies and a shared, hot 
intragroup medium. 



2. Sample Definition, Observations, and Data Reduction 

Objects for study were selected on the basis of data available from the archive at the Chandra 
X-ray Center (CXC) for CGs that are part of a sample designed to study star formation and galaxy 
evolution in the compact group environment (see Johnson et al.|20"07 for more information regarding 



the original sample; HCG 92 was subsequently added based on the availability of high-resolution, 
deep imaging data). To this end, the sample consists of groups at small enough distances (z < 0.015, 
excepting HCG 92) to allow for high spatial resolution imaging over a wide range of wavelengths. 
In addition to our new Chandra observations of HCGs 7, 22, 31, and 59 (Pis: Garmire [7, 22] and 
Gallagher [31, 59]), we selected archival data that covered the entirety of each group, rather than 
only individual group members, to search for hot, ionized gas in the intragroup medium, bringing 
the total number in our study of diffuse X-ray emission to nine groups. HCGs 2, 19, 48, and 61, 



which were included in the Johnson et al. (2007) sample, were omitted from our study due to the 
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lack of suitable Chandra imaging data. 

All of the groups were observed at the ACIS-S aimpoint with the exception of HCG 90, which 
was observed at the ACIS-I aimpoint due to its large angular extent (the diffuse emission of HCG 90 
falls mostly on the 10 CCD, however some is present along the edges of the other I array detectors 
as well). The individual ACIS CCDs have a field of view of 8f 3 x 8(3, which is comparable to 



the extent of the diffuse emission in compact groups found in Ponman et al. (1996) and Mulchaey 



et al.| Q2003| ) (noting the exceptions of HCGs 62 and 90). All data were taken in VFAINT mode 



except HCG 16, HCG 62 (obsID 921), and HCG 92 (obsID 789), which were in FAINT mode. 
Due to the design of the High Resolution Mirror Assembly, observations with the Chandra ACIS 
camera have superior spatial resolution (~0.5" FWHM) compared to other instruments such as the 
ROSAT PSPC (~25 /; ) and XMM EPIC (~6"), allowing for more robust removal of point sources 
that contaminate the diffuse emission and can ead to incorrect estimations of hot gas properties. 
Tables [3] and [T] list the properties and observation information, respectively, of the CGs included 
our sample. 

The data were reprocessed beginning with the level 1 events file. We used the CIAO tool 
acis_reprocess_events with standard event filtering and VFAINT background cleaning (when 
appropriate) to produce the level 2 events file; however, we omitted the pixel randomization step 
in the event file processing to prevent degradation of the spatial resolution. The pixel randomiza- 
tion introduces a G/.'5 random resampling of the event positions on the detector. This procedure 
is normally performed to mitigate the effects of aliasing in observations < 2 ks, however the ob- 
servations included in our sample are substantially longer than this limit, and therefore the pixel 
randomization is not required. The pixel randomization has been shown to decrease the spatial 
resolution of the observations by ~12% (e.g., Chartas et al.||2002 ). High spatial resolution is useful 



for proper removal of point sources that may be embedded within or projected onto any diffuse 
emission present in the groups. 

Because the diffuse intragroup emission is typically weak and there is no robust method to 
determine the full radial extent of it from the data, for most targets we did not perform a spectral 
extraction in an annulus surrounding the targets to determine the soft X-ray background level. 
Instead, we subtracted the instrumental background using the stowed ACIS background files ap- 
propriate to each observation (as in White et al. 2003). The stowed ACIS background files^were 



obtained by sliding ACIS out of view of the sky and away from the external calibration source to 
avoid spectral-line contamination. All stowed background observations were taken after 2002 in 
VFAINT mode with gain corrections applied. Following procedures outlined in the ACIS Extract 
manual ( Broos et al.||2010 ), we removed the "Clean55" bit from the event lists for subtraction from 



data taken in FAINT mode and renormalized the background files to the data using the particle 
background in the range 9-12 keV. The stowed background files include emission from the high 



1 This brief description of the stowed ACIS background is based on the work of Maxim Markevitch; further details 
are at http://cxc.harvard.edu/contrib/maxim/stowed/ 



Table 1. Chandra Observation Information 



HCG Chandra 


Date 


Extraction Area a 


Total Counts' 3 


Exp. 


Bkg Counts 


Bkg Exp. c 


obsID 




(sq. arcmin) 


(0.7-3.0 keV) 


(ks) 


(0.7-3.0 keV) 


(ks) 



7 d 


8171 d 


2007-09-13 


48.2 


2419 


19.2 


57670 


900 




9588 d 


2007-09-16 






16.7 






16 


923 


2000-11-16 


26.9 


1270 


12.6 


6209 


110 


22° 


8172 


2006-11-23 


44.7 


2078 


31.8 


27973 


450 


31° 


9405 


2007-11-15 


8.70 


636 


35.6 


5433 


450 


42 


3215 


2002-03-26 


36.5 


7600 


31.7 


15882 


450 


59° 


9406 


2008-04-12 


24.7 


1457 


38.4 


15373 


450 


62 


921 


2000-01-25 


7.54 


31053 


48.5 


2604 


145 




10462 d 


2009-03-02 


7.54 


55302 


67.1 


247702 


900 




10874 d 


2009-03-02 






51.4 






90 


905 


2000-07-02 


68.3 


24821 


49.5 


166640 


450 


92 


789 


2000-07-09 


19.3 


2545 


19.7 


4663 


110 




7924 


2007-08-17 


19.3 


10656 


93.2 


27337 


450 



Notes. 

a Total area over which the events were extracted 
b Source + Background in science observation 
c Exposure time of the ACIS stowed background files 
d Count information listed for merged spectra 



c Counts evaluated over 0.7-1.7 keV (see E 3.2 ( 
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energy particle background (E > 2 keV), but not the soft X-ray background (SXRB; E < 2 keV), 
therefore in groups where we used the stowed background files we expect contamination from the 
SXRB at energies <2 keV. Table [T] lists total counts extracted from and exposure times of the 
stowed background data for each group. 

We used the latest version of ACIS Extract to create exposure maps, model point source PSFs 
and excise them from the observations, and extract spectral information (see ^3]). The energy used 
to create the exposure maps was the median event energy over the range 0.5-3.0 keV, chosen to 
coincide with the peak of any soft diffuse emission in the data. 



3. Spectral Extraction and Modeling 
3.1. Point Source Detection and Removal 



We used the CIAO Mexican-hat wavelet detection routine wavdetect (Freeman et al. 2002) 
to search for point sources in the field of view down to a source-significance threshold of 10 -5 . 
Detections were examined by eye to remove spurious sources (e.g., those with an axis length < 
1 pixel). We then used ACIS Extract to model the PSF of the sources with MARX prior to 
excising the point sources from the events files. The excised regions correspond to 1.1 times the 
99% encircled energy radius. For extraction of spectra, we did not interpolate over the holes 
created by excising the point sources to avoid making statistical assumptions concerning the gas. 
The area of diffuse emission on the sky lost due to point-source removal was typically less than 
1% of the total extraction region. We defer the examination of the point sources in our sample to 
Tzanavaris et al. {in prep). 



3.2. Extraction of the Diffuse Emission 



Extraction regions were centered on the apparent centers of the CG galaxy distributions with 
shapes and sizes chosen to best cover all of the main group members on the S3 chip (10 for HCG 90), 
including evident diffuse emission in the level 2 events files. The extraction regions were either 
circular or elliptical in shape, except in the case of HCG 16, which fills approximately half of the 
S3 CCD. In this group, the region was rectangular and placed at an angle to include the main 
galaxies and the known distribution of H I gas. The area of each extraction region can be found in 
Table [TJ while additional information is located in Appendix [Aj Note that in the case of HCG 22, 
the extraction region includes a background pair of galaxies to the southeast of the group center. As 
this group was classified as a non-detection, we set an upper limit on the diffuse X-ray luminosity 
using the method outlined below. Therefore, the inclusion of this pair in the extraction region has 
negligible impact on the results. 



Prior to model fitting in XSPEC (Arnaud 



1996), the extracted spectra were rebinned such 
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that each pulse height amplitude (PHA) bin contained a minimum of 20 counts. This ensured 
statistically valid results when fitting with x 2 statistics. In cases where the ACIS instrumental 
response did not change significantly between observations (i.e., multiple obsIDs within an observing 
cycle), we merged the on-source spectra to reduce the relative error in each PHA bin. When merging 
the spectra was not possible due to the changes in the instrument response, we simultaneously fit 
the model to each spectrum with the temperature and metal abundances linked between different 
obsIDs; however, the model normalization was allowed to vary freely. 

A detection was defined to have S/N > 3. The noise of the spectra was determined to be 

1/2 

a = [SB + (A s t s / A^t^B] 1 , where SB is the total counts in the source before background sub- 
traction, B is the number of counts in the background, A is the area of the extraction region, t 
is the integration time, and the subscripts s and b represent the science and background observa- 
tions, respectively. When using the stowed background, the background counts were scaled to the 
9-12 keV count rates in the science observations prior to estimating the noise. This scaling has the 
effect of normalizing the shape of the stowed background spectrum to that of the source in a range 



of energies dominated by the particle background Hickox & Markevitch (2006). When no diffuse 
emission was detected, we put an upper limit on the luminosity of the intragroup medium by fixing 
the temperature and metal abundance of the plasma to reasonable values of 0.6 keV and 0.5 Z Q , 
respectively. We then adjusted the normalization to match the sum of the observed count rate and 
the la noise estimate. 

In all cases when diffuse X-ray emission was not detected, the energy range over which the 
count rates were evaluated had to be restricted to 0.7-1.7 keV due to oversubtraction of the stowed 
background between 1.7 and 2 keV. Bolometric X-ray luminosities were then computed in the same 
manner as for groups with detected emission (see below). 



3.3. Spectral Model Fitting 



We fit a combination of optically thin plasma and Galactic photoelectric absorption models for 
each of the groups. The best fitting models are presented in Table [2] The first model component, 



an optically thin plasma, was modeled using the MEKAL model (Mewe et al. 1985 1986 Kaastra 



1992, 1993 Liedahl et al.||1995 Kaastra Sz Liedahl|1995 ) with the adopted ionization balance taken 



from Arnaud & Rothenflug (1985) and Arnaud Sz Raymond (1992). We fixed the hydrogen number 
density at a reasonable value of 1 cm -3 and allowed the remaining parameters to vary freely. While 
the metallicity parameter was allowed to vary, these values are poorly constrained by the available 
data, and we only report them here to describe the best-fitting models to the observations. We 
calculated the model at all temperatures rather than interpolating it from a pre-calculated table 
while fitting the data. Other than the plasma temperature, we also report the normalization used in 
calculating the X-ray luminosity. We note that because the gas is often associated with individual 



galaxies rather than a single distribution permeating the intragroup medium (see £4.1), it is likely 
multi-temperature; however, we find that the single temperature plasma model fits the extracted 
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spectra well, except in the case of HCG 62 for which a two-temperature plasma results in a better 
fit. Note that HCG 62 has the most counts in the background-subtracted data compared to any 
other observation in this study; specifically, a factor of ~7 more counts in the merged obsIDs 10462 
and 10874 compared to HCG 90, which has the next most counts. For other targets, our ability to 
isolate the multiple spectral components that likely make up the diffuse X-ray emission is limited 
by the number of counts; therefore, it is probable that we cannot detect spectral complexity in our 
data (e.g., multi-temperature plasmas) given the limited numbers of counts. 

Photoelectric absorption was modeled using the Tuebingen-Boulder ISM absorption model 



from Wilms et al. (2000). The only model parameter, the H I column density along the line of 



sight, was fixed at the value determined using the HEASoft tool nH to compute the weighted mean 
of H I in a cone centered on the source and with a radius of 1°. Following the recommendation of 
the nH manuaj^J we use the H I values from Kalberla et al. ( 2005 ) . 



Abundances and depletion values for the model (used in both the absorption and emission 



components) were taken from Lodders (2003) (the most recent abundances available in XSPEC) 



rather than the default from Anders & Grevesse (1989). Both abundance tables use Solar photo- 



spheric line and CI chondrite analyses to determine the relative amounts of each element. We found 



that the values from Lodders (2003) consistently performed better at fitting emission-line features 



in the spectra, particularly in observations with significant numbers of counts (e.g., HCG 62). 

Because Chandra only has significant response over the energy range 0.3-8 keV, we computed 
the bolometric X-ray luminosity Lx using a dummy response over the energy range 0.01-100 keV 



logarithmically divided into 5,000 energy bins. We used 3K CMB adjusted velocities (Fixsen et al 



1996) from the NASA Extragalactic Database (NED), which adjusts the velocities for the observed 
dipole anisotropy in the CMB, to determine the distances to the sources for luminosity calculations. 
Parameters of interest, specifically the temperature and luminosity of the plasma (and by necessity, 
the model normalization), are reported with 90% confidence error bars. For simultaneously fit 
spectra, we report the average luminosity (weighted by the number of counts) determined from 
the model fits. Note that in the case of HCG 59, the peak in energy of the X-ray emission is very 
poorly constrained; therefore we fix the temperature to the best-fitting value prior to determining 
the error in the model normalization. Furthermore, we do not consider the best-fitting value to 
be representative of the real temperature and simply report it as < 1 keV, however we do use the 
temperature result of the model fit for qualitative purposes in the figures below. We present the 
temperatures and luminosities derived from the fitted models, as well as goodness of fit estimations, 
in Table EJ 



http : //heasarc .nasa.gov/Tools/w3nh_help . html#compari son 
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Table 2. Best Fit Spectral Model Parameters 



Tjppa 

1 1 V_ V I 


W T 


K 1 


z/z 


Amekal 


!„„ ( T ,,\ 
iOgio^Xj 


,,2 /pi n TT 




(10 20 cm- 2 ) 


(keV) 




(10~ 4 cm- 5 ) 


(erg s -1 ) 




7 b 


1.97 


0.6 


0.5 


< 0.29 


< 40.35 




16 


2.56 


0.65 + 0.06 


0.16 


3 Q0 +2 ' 20 


41 10+° 19 


40.43/47 


22 b 


4.26 


0.6 


0.5 


< 0.21 


< 39.84 




31 


5.70 




0.05 


1 16+ 3 ' 48 

1 ' 1D -0.65 


40-54±g:|g 


31.43/40 


42 


4.11 


0.72 + 0.01 


0.55 


8-38±i;?l 


41.82±8:?g 


196.48/114 


59 


2.64 


29+ - 41 


17.53 


0.01 ± 0.004 


40.08t° : l 2 


55.84/61 


62 (921) cool 


3.32 


0.71 + 0.01 


16.68 


n c .n+ 1 - 23 


41.88;?; 54 


319.10/281 


62 (921) hot 


3.32 


1.20 + 0.04 


0.48 


1-461^ 


42.00 ± 0.04 


319.10/281 


62 (921) total 










42 23+ - 34 

4lZ.ZO_o 26 




62 (10462+10874) cool 


3.32 


0.71 + 0.01 


16.68 




41.84;?' 54 


319.10/281 


62 (10462+10874) hot 


3.32 


1.20 + 0.04 


0.48 




42 ni+° : ° 4 


319.10/281 


62 (10462+10874) total 










42 23 +0 ' 31 




62 average total 










/10 79+0-33 
^Z. (Z_ 26 




90 BCD 


2.02 


0.66 + 0.03 


0.54 


2 66+ 159 


40.79;g 2 g 


121.75/119 


92 (789) 


6.16 


0.63 + 0.02 


0.14 


6 25+ 1 ' 08 

°- zo -0.99 


+0 :o7 

' ' -0.08 


357.33/220 


92 (7924) 


6.16 


0.63 + 0.02 


0.14 


5.8i±l:?l 


41 7Q+0.09 
J -0.10 


357.33/220 


92 average 










41.76+0.07 





Notes. 

a For groups with spectra that were fitted simultaneously rather than merged, we list the corresponding 
obsIDs in parentheses 

b No detection, model temperature and abundance fixed at 0.6 keV and 0.5 Z Q , respectively 
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Table 4. Comparison of Chandra and R OSA T Measurements 



HCG 


Tx,P96 


T x 


AT X 


-fx,P96 a Fx 3 " 


A Fx 




(keV) 


(keV) 


(keV) 


(10~ 12 erg s _1 cm -2 ) (10 -12 erg s _1 cm -2 ) 


(10~ 12 erg s" 1 cm" 2 ) 



16 


0.30 ±0.05 


0.65 ±0.06 


-0.35 ±0.08 


42 


0.82 ±0.03 


0.72 ±0.01 


0.10 ±0.03 


62 b 


0.96 ±0.04 


0.99 ±0.04 


-0.03 ±0.06 


90 


0.68 ±0.12 


0.66 ±0.03 


0.02 ±0.12 


92 


0.75 ±0.08 


0.63 ±0.02 


0.12 ±0.08 



7.97±i 
23.71 
169.48j 



18 
03 
1-1.12 
-1.07 
12.12 
11.31 

11.34i£S 



9.04 



+0.87 
0.80 



-1.34 
-0.99 
H4.38 
-4.34 

158.30™ 

4 47+ 2 ' 67 
*- 4 '-2.53 
+ 1.20 

1.09 



1. 

20. 74^ 



6.95 



5.99Z 



11.18 



1.78 
1.43 

2 Q7+ 4 - 52 
' -4.47 
182.59 



-71.' 



2.09li;f 5 



Notes. 

P96 = |Ponman eTal] ( [1996 1 

a Converted from the bolometric luminosity and adjusted for differences in cosmology 

b The values listed for this study are the flux-weighted average temperature and total flux of the two plasma components in 
HCG 62 
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Results and Discussion 



We identify hot diffuse gas in seven of nine groups (HCGs 16, 31, 42, 59, 62, 90, and 92) 
in our X-ray sample. Of the detected groups, diffuse emission in HCG 31 has not been reported 
in previous studies. Both of the two groups without detections, HCGs ij^] and 22, have low to 
negligible star formation without much evidence of strong tidal interactions in the past few Gyr 



(e.g., Konstantopoulos et al. 2010). In the detected groups, the temperatures are all fairly similar 



(0.6-0.72 keV not including the hot component of HCG 62 or the anomalously low temperature of 
HCG 59, which is very poorly constrained), while the range in X-ray luminosities spans 10 40 - 
10 42 18 erg s" 1 . The morphology of the hot gas ranges from isolated around individual group 
members to common X-ray halos (e.g., HCGs 16 and 62, respectively), as well as gas bridges 



connecting galaxies (HCGs 59 and 90; see £4.1). In Figures [lj|7| we plot all values derived from 
model fits for multiple observations that could not be merged prior to fitting. 

Table [3] lists relevant information about each of the groups in our sample. Included in the table 
are the group redshift and 3K CMB velocity, the calculated velocity dispersion, number of main 
group members, the number of galaxies used in the velocity dispersion calculation, the number of 
E/S0 type galaxies that are not considered dwarf members, the specific star formation rate from 



Tzanavaris et al. (2010), the total H I mass from Verdes-Montenegro et al. (2001) and Borthakur 



et al. (2010), the dynamical mass determined from the group velocity dispersion and mean two 



galaxy separation taken from Hickson et al. (1992), and the group H I evolutionary type defined by 



Johnson et aL] fl2007| ) (see £ |4.3| for a full description) . We will continually refer back to these data 
throughout the figures and discussion that follow. 

The nine groups in our sample were also observed with ROSAT and the results presented in 
Ponman et al. (1996). We compare our temperatures and luminosities derived from fitting the 



extracted spectra from the Chandra data with those also detected by Ponman et al. ( 1996[ ) (i.e., 
HCGs 16, 42, 62, 90 and 92) in Table[4} Rather than directly compare the luminosities, we calculate 
the fluxes of the sources, which removes the dependence of the assumed cosmology. The redshifts 
used in the flux calculation for the ROSAT data are the same as those listed in the redshift column 
of Table [3} Note that [Ponman et ah] ( fl996| ) give errors in the luminosity and temperature at the la 
level. For the temperature in HCG 62, we use the luminosity-weighted average of the two plasma 
temperatures, while the fluxes in HCGs 62 and 92 are based on the average of the luminosities 
weighted by total counts in the different Chandra obsIDs, therefore the percent errors on these 
values seem quite large because of the compounding of errors from multiple measurements. In 
all cases, the luminosity of the X-ray emission from the ROSAT data is found to be brighter in 



3 A deeper 49 ks XMM observation (PI: Belsole) of this group exists, but no results have been published at the 
time of this writing. A brief inspection of the data shows several dozen point sources and large-scale diffuse emission 
across the field of view of the EPIC camera (~5 times the angular size of HCG 7). The sensitivity of XMM to cool, 
extended gas complicates the interpretation of this emission as either related to the group or part of the SXRB. 
Results from this and other XMM data will be included in subsequent papers. 
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the Ponman et al. (1996) study. When comparing the temperatures we determined against those 
from Ponman et al. (1996), we find that HCGs 62 and 90 are consistent between the two studies, 
while HCG 16 is significantly hotter, and HCGs 42 and 92 are significantly cooler than previously 
reported. The mean temperature of the groups detected in both studies is 0.70 and 0.72 keV with 
standard deviations of 0.25 and 0.14 keV in Ponman et al. (1996) and this study, respectively. The 
temperature discrepancies are likely due to the difference in the PSF of ROSAT PSPC instrument 
compared to the Chandra ACIS camera. Specifically, the detection and subtraction of point sources 
in the diffuse emission is much more robust using Chandra data, and point sources contaminating 
the extracted spectrum would alter the peak of the emission leading to an incorrect estimation 
of the hot gas temperature. We also note that Ponman et al. (1996) attempted to subtract the 
diffuse X-ray contribution of the individual group members and interpolate over the 'holes'. This, 
combined with distinctions in the extraction regions and the responses of the ROSAT PSPC and 
Chandra ACIS instruments, may explain the differences in the results. 



4.1. Gas Distribution and Morphology 

Figure [T] shows the smoothed diffuse X-ray emission in the CGs included in our study. To 
construct the X-ray maps, we smoothed the level 2 events files in the range of 0.5-2 keV using 
the adaptive smoothing csmooth algorithm. Pixels with S/N > 3 above the stowed background 
were smoothed with a Gaussian kernel with smoothing scales between 2 and 10 pixels. From the 
figures, we note that the hot gas in compact groups exhibits varied morphologies including small 
halos around one or several galaxies, plumes centered on particular group members, tidal bridges, 
and large common halos encompassing most of the group. 

The maps were qualitatively assessed by eye to categorize the observed X-ray emission as either 
associated with the environment or the individual group members. Of the nine CGs included in 
our study, we do not detect diffuse emission from HCGs 7 and 22. From those groups detected with 
Chandra, only HCG 62 has emission that permeates the IGM, while HCG 42 has a bright X-ray 
halo centered on the brightest group galaxy (i.e., 42A). We consider these two groups indicative 
of systems with X-ray emission similar in morphology to clusters. The other groups, HCGs 16, 
31, 59, 90, and 92 (excepting the shock front in the case of HCG 92), are those in which the 



hot gas is associated with the individual galaxies. Mulchaey (2000) noted that previous X-ray 
telescopes did not have the requisite spatial resolution to separate the intragroup gas from the 
galaxy-linked emission; however, the resolution of Chandra is well suited to this task. Futhermore, 
we join Tamburri et al. (2012), who presented a Chandra study of HCG 79 (also known as "Seyfert's 
Sextet"), in speculating that perhaps many compact groups exhibit galaxy-linked X-ray emission 
rather than a hot intragroup medium associated with the environment. Throughout the remainder 
of the paper, unless explicitly stated, we discuss the galaxy- and environment-linked diffuse X- 
ray emission together as they fit into a larger picture of group evolution, i.e., from dynamically 
unevolved systems with X-ray emission confined to the galaxies to more evolved systems with a 
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single X-ray halo. We stress that it is difficult to disentangle the galaxy- and group-linked emission 
in the data as this distinction is not always clear, as in the emission surrounding HCG 42A. 




40.0 30.0 0:39:20.0 10.0 

Right ascension 



Fig. 1. — Smoothed 0.5-2.0 keV X-ray contours of groups in our sample. The event files were 
smoothed using the CIAO task csmooth with the minimum significance set to 3cr above the stowed 
background scaled by the 9-12 keV count rate. The minimum and maximum smoothing scales used 
were 2 and 10 pixels, respectively. In all images, the dashed blue line and red crosses correspond 
to the extraction region used (see Appendix [AJ and the locations of the X-ray point sources that 
were excised from the analysis (but retained in the smoothed X-ray data), respectively. The optical 
images are DSS POSS2 survey red filter data. Contour levels in counts follow in the captions for 
each group. HCG 7: 0.1, 0.2, 0.3. 
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Fig. 1 — Continued. HCG 16: 0.1, 0.5, 1, 3, 5, 10. HCG 22: 0.1, 0.2, 0.3. 
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Fig. 1.— Continued. HCG 31: 0.2, 0.5, 1, 3, 5. HCG 42: 0.5, 1, 3, 5, 10. 
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Fig. 1.— Continued. HCG 59: 0.2, 0.5, 1, 3, 5, 10. HCG 62: 3, 6, 10, 15, 30, 50. 
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Fig. 1.— Continued. HCG 90: . HCG 92: 0.5, 1, 2, 5, 10. 
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HCG 16 




» 



Fig. 2. — Three color image of HCG 16. Red and green represent Ha and i?-band data from CTIO 



(Meurer et al. 2006), respectively, while blue is the smoothed soft (0.5-2.0 keV) X-ray emission 
from Chandra. North is up and East is to the left in the image. The scale bar indicates 1'. The pink 
color in the nucleus of galaxy A and surrounding galaxies C and D is caused by the overlapping 
of Ha and X-ray emission. Galaxy D exhibits a biconical gas distribution centered on its nucleus, 
while ionized gas surrounds galaxy C with an elongation along the North-to-South axis. 



In the specific case of HCG 16, the group contains two members (C and D) with hot gas 
plumes that are coincident with extended Ha emission (see Figure . Both galaxies show clear 



signs of recent interactions (e.g., disturbed velocity fields; gas/stellar disk misalignments; Mendes 



de Oliveira et al. 


1998 


2001 


) and are strong 



respectively, and SSFRs of 



(2 - 3) x 10~ 10 yr" 1 ( |Ribeiro et aL||l996t |Tzanavaris et al.||2010[ ). |Rich et al.| ( |2010[ ) used integral 
field spectroscopy to demonstrate that the motions of the gas and the optical emission-line profiles 
of 16D are indicative of an M82-like superwind due to the intense starburst within the galaxy, 
and that 16D and M82 share some similar properties (e.g., metallicity, Ljr, Ln a ). Based on its 
high star formation rate and extended emission-line gas, 16C may represent a second M82-like 
starburst in the HCG 16 group, possibly triggered by a recent (< 1 Gyr) tidal interaction with 
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16A similar to the encounter between M81 and M82 ~220 Myr ago ( Gottesman fc Weliachew|1977 



Konstantopoulos et al. 2009). Furthermore, Jeltema et al. (2008) found evidence for a weak soft 



X-ray bridge connecting galaxies 16A and B indicating a recent tidal interaction between these 
group members as well. 

Further examining the morphology of the hot gas and its relation to the galaxies in the compact 
groups, we find that the gas in low dynamical mass (< 10 12 ' 5 M ) groups with low velocity disper- 
sions (< 250 km s _1 ) is concentrated around the individual group members, while gas in relatively 
higher mass systems with larger velocity dispersions begins to resemble a common envelope. 



4.2. The Lx — T and Lx — cr Relations Re-visited 



Previous work has examined the relationships between the X-ray luminosity and both the 



plasma temperature and the cluster/group velocity dispersion in systems of galaxies (e.g., Ponman 
et~aI1[T996] |Mulchaey fc Zabludoffj [T9981 |Wu et al.|[i~999t |Helsdon fc Ponman| [2000} |Osmond fc 
Ponman]|2004 Mulchaey et al.||2003 ). The cluster data have been found to show very little scatter 
in Lx — T and Lx — cr space. However groups, with fainter X-ray luminosities, have been found 
to exhibit a larger spread potentially due in part to uncertainties in the measurements, or because 



they are not virialized systems (see below). We merge the cluster data from Wu et al. (1999) 



and Zhang et al. (2011) for comparison to the CGs. Prior to merging, we first adjusted the Wu 



et al. (1999) data to our assumed cosmology, while Zhang et al. (2011) used the same cosmology 
as that assumed in this paper. When comparing the cluster data to the compact groups in Lx — T 
space, we only include the 176 clusters with uncertainties in both Lx and temperature. This same 
criterion, with respect to velocity dispersion rather than temperature, is applied to the cluster data 
when comparing the groups and clusters in Lx — cr space, resulting in 142 clusters in this sample. 



We note that the Zhang et al. (2011) clusters are measured to rsoo (i-e., the radius at which the 
average density falls to 500 times the critical density at that redshift), and the extraction regions 



used in this work vary between 26% and 69% of rsoo- It is unclear to what radius the Wu et al 



(1999) clusters were measured, however the authors state that they used a /3 model to correct all 



of their data to the same fraction of the virial radius. 



Figure [3] shows the Lx — T relationship for the groups and the clusters from Wu et al. ( 1999 ) 



and Zhang et al. (2011 ). Using linear regressions described by Akritas & Bershady ( 1996 ), we fit the 



clusters in Lx — T space not including the compact groups. For this fit, we use the orthogonal dis- 
tance regression (ODR) fitting method, which finds \og w {L x ) = (42.2±0.2)+(3.33±0.20) log 10 (T), 
where Lx and T are in units of erg s _1 and keV, respectively. Examining Figure [3J we find that 
HCGs 31, 42, 59, and 92 agree with the Lx — T cluster relation within the errors. We note that 
HCG 62 lies slightly above the the cluster relation (i.e., it is brighter for its temperature). HCGs 7, 
16, 22, and 90 lie below the Lx — T relation from the clusters. 
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-0.5 0.0 0.5 1.0 

Group X — ray Temperature log, (T) [keV] 



Fig. 3. — The Lx — T diagram for compact groups (black filled circles). The solid line indicates the 
best fit to the cluster data (gray crosses) from Wu et al. (1999) and Zhang et al. (2011) using an 
orthogonal distance regression (ODR) fit to 176 clusters, while the dashed lines represent the errors 
on the fit. We plot the two plasmas in HCG 62 as one data point using the luminosity- weighted 
average of the temperatures from Table [4] and the total X-ray luminosity. The temperature of 
HCG 59 is very uncertain due to the poorly defined peak in the X-ray spectrum, and we conserva- 
tively estimate it to be <1 keV, however we include the best-fitting temperature from the MEKAL 
model here in the Lx — T diagram for qualitative purposes. The best fitting cluster relation agrees 
with HCGs 31, 42, 59, and 92 within the errors, while HCG 62 agrees lies within the scatter in the 
cluster data. We do note most of the compact groups lie systematically below the cluster fit. 
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Fig. 4. — The Lx — <? diagram for compact groups. Due to the inaccuracy in the velocity dispersions 
for systems with a small number of components, we plot filled black circles to indicate groups 
for which 6 or more galaxies were used to compute the velocity dispersion, while open triangles 
represent groups below this threshold. The solid line indicates the best fit to the clusters using a 



linear-least squares fit to the 142 clusters taken from Wu et al. ( 1999 ) and Zhang et al. (2011 ) (gray 
crosses), while the dashed lines represent the errors on the fit. The best fitting cluster relation 
agrees with HCGs 16, 31, and 62, while HCGs 7, 42, 59, 90, and 92 lie below the fit. We note that 
the agreement of HCGs 16 and 31 with the cluster fit is likely coincidental because the emission in 
these systems is largely due to star formation rather than virialization, therefore they should not 
necessarily agree with the fit to the cluster data. 
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Figure [4] shows the Lx — cr relation for the groups in our sample. Velocity dispersions were 
calculated using the most accurate velocities available from NED for the CG members including 



additional known group members primarily from the large spectroscopic surveys of de Carvalho 
etll| ( p97l ), |Zabludoff Mulchaey| ( fl998| ), and |Zabludoff fc Mulchaey| pOOO] ) (see Table [| for 
full velocity references). Velocity errors for individual group members were typically < 30 km s _1 . 



Again, we compare our compact group data with the cluster data taken from Wu et al. ( 1999 ) and 



Zhang et al. (2011) using a linear least-squares fit. We chose this type of fit because the ODR 
method used in the Lx — T diagram showed strong systematics in the residuals of the cluster data. 
The resulting best fit is log 10 (Lx) = (33.7 ± 0.6) + (3.70 ± 0.21) log 10 (a), where Lx and a are in 
units of erg s _1 and km s _1 . We find that two of our CGs (HCGs 16, 31) agree with the Lx — cr 
relation from the galaxy clusters within errors, while HCG 62 is similar to the clusters within the 
scatter. In the cases of HCGs 16 and 31, this agreement is likely coincidental because the Lx — cr 
relation is predicted from the Virial theorem to be Lx oc <r 4 , and the bulk of the hot gas in these 
systems is clearly due to star formation. We do note, however, that for groups in which dynamical 
processes are increasingly important in heating the gas (i.e., HCGs 42, 59, 62, 90, and 92), there 
does appear to be a monotonic increase in the X-ray luminosity with increasing velocity dispersion, 
albeit with some scatter. 



Ponman et al. (1996) hypothesized that higher temperatures for a given Lx than would be 
expected from the cluster Lx — T relation for compact groups could be a result of the injection of 



energy into the IGM from galactic winds, while Ponman et al. (1999) suggested that the observed 



deviation of the galaxy group X-ray luminosities from the established cluster scaling relations could 
be explained by preheating of the IGM by supernovae. In the preheating model, heating of the gas 
in the IGM occurs early in the lifetimes of groups and similarly steepens the Lx — T relation for 
these systems. Despite groups such as HCG 16 where star formation may play a more prominant 
role, the existence of X-ray brighter groups (i.e., HCGs 42) below the Lx — T scaling relation for 
clusters indicates that galaxy winds alone are unlikely to explain the observed deviation. However, 
the observed dichotomy in X-ray gas morphology (i.e., galaxy- vs. environment-linked emission) 
may indicate that the low- mass, low velocity dispersion systems are dynamically unevolved and 
have shallow potential wells that are unable to heat any neutral gas that has been liberated from 
the group members. In these systems, the role of individual galaxies may be more important 
in heating gas through local (e.g., star formation, superwinds, and accretion) rather than global 
processes (e.g., virialization) . As systems accrete more mass or as the system relaxes, the potential 
well will deepen and the contribution of individual group members to the diffuse X-ray emission 
should lessen. 

After attempting to subtract the X-ray emission associated with the group members in HCGs, 
Ponman et al. (1996) found that the Lx — cr relation was flattened compared to the clusters (i.e., 



the rate of change in Lx as a function of a was slower for groups than for clusters), similar to the 
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delP Antonio et al. (1994) Einstein study of rich groups^} as well as a subsequent study of CGs by 
Helsdon &; Ponman (2000). The morphology of the hot gas in CGs may provide an explanation 
for the groups that disagree with the Lx — T and Lx — cr relations derived from galaxy clusters. 
As discussed in §4.1[ the hot IGM that permeates HCG 62 and the X-ray bright halo centered on 
HCG 42A have morphologies that are qualitatively similar to the hot gas observed in the cluster 
environment, i.e., where the gas has been heated by virialization. In the remaining groups, the only 
diffuse X-ray emission we detect is associated with the individual galaxies rather than the group 
environment itself. This suggests that the gas is not in hydrostatic equilibrium within the group, 
and therefore the temperature and luminosity of the total X-ray gas in the system does not trace 
the group potential as it does in the galaxy clusters. 

For completeness, we compare our sample of compact groups to studies of the X-ray scaling 
relations in both normal (i.e., those that are not compact) and fossil groups. With respect to the 
Lx — T relationship, Khosroshahi et al. (2007) and Harrison et al. (2012) noted that fossil groups 



match well with both clusters and normal groups. Only three of the compact groups presented in 
our study (HCGs 42, 62, and 92) agree well with the cluster scaling relation within the scatter, 
though HCG 92 is dominated by emission from a strong shock. The remaining CGs at T < 0.65 keV 
fall below the best fit to the cluster data. In Lx — cr space, Rines & Diaferio (2010) showed that 



normal groups did not deviate from the relation for clusters, while Khosroshahi et al. (2007) found 



that fossil groups are consistent with the best fit to the cluster data, though more X-ray luminous 
than the normal groups in their sample. Our sample of compact groups fall systematically below the 
observed Lx — cr relation for clusters with the exceptions of HCGs 16 and 31, which are dominated 
by vigorous star formation. This may be in agreement with the interpretation of merging systems 



being X-ray underluminous for their velocity dispersions (e.g., Rasmussen et al. 2006 Popesso 



et al. 2007). We caution the reader that in most cases the velocity dispersions of the groups 



are measured from very few galaxies, therefore they likely do not accurately represent the three 
dimensional dispersions that are assumed by the Lx — cr relation. 

Based on the discrepancies between the compact groups and the X-ray cluster scaling relations, 
we postulate that systems similar to the low-mass, X-ray faintest groups in our study should not be 
considered analogs to clusters, with the possible exception of HCG 62, which lies close to the cluster 
data in both Lx — T and Lx — cr within the scatter. It is important to note that this dissimilarity 
between the low-mass groups and clusters does not preclude these systems from becoming more 
cluster-like if they somehow become similar to the more massive, rich groups in our sample (e.g., 
by accreting additional members and continuing strong interactions to liberate gas into the IGM). 
However, it is unlikely that rich groups today formed from poor groups like those observed in 
the current epoch. Therefore the exact mechanism whereby the poor groups in our sample could 
become more cluster-like remains unknown. We do note, however, that there could be extended, 



4 While dell' Antonio et al. ( 1994 \ did not subtract the galaxy-linked emission from their data, they did avoid using 
data in instances where the hot gas was clearly associated with only galaxies and not the group environment. 
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faint emission that is undetectable in the available data (e.g., the X-ray IGM of HCG 16; Belsole 



et al. 2003). Deeper observations may reveal cooler, lower luminosity gas associated with entire 



groups rather than the galaxies, which could shed light on the relationship between poor groups 
and more cluster- like systems. 



4.3. Comparison of the X-ray Data with H I Gas 



Previous work has shown that relatively H I rich compact groups contain galaxies that ex- 
hibit mid-IR colors dominated by star formation with correspondingly high SSFRs measured from 



UV+24 fim fluxes, while groups deficient in H I have more quiescent colors and low SSFRs (Johnson 



et al.|2007 


Walker et al.|2010 


2012 


Tzanavaris et al.|2010 



of compact groups in which the H I gas is processed either through star formation or ionization by 



the group potential (the velocity dispersions imply virial temperatures of ~0.08-0.3 keV; Verdes 



Montenegro et al. 2001 Johnson et al. 2007). Therefore, we expect that the most H I abundant 
groups should have very weak diffuse X-ray emission, while the H I poor groups should have the 
brightest X-ray luminosities. Furthermore, the morphology of the H I gas may dictate how it is pro- 
cessed, i.e., through star formation or by virialization in the IGM in CGs with neutral gas confined 
to the galaxies or stripped into the IGM, respectively (cf. Konstantopoulos et al.|[2010 ). Through- 
out this section, we use the H I abundance type notation from Johnson et al. (2007). Types I, 



II, and III indicate decreasing H I abundance relative to the group dynamical mass, respectively. 



Johnson et al. (2007) quantitatively define these H I mass types as (I) log(Mn i)/log(Md yn ) > 0.9, 
(II) 0.9 > log(M H i)/log(M dyn ) > 0.8, and (III) log(M H i)/log(M dyn ) < 0.8. The H I evolutionary 
types for each group are listed in the last column of Table |3j 

We calculate the H I to dynamical-mass ratios for the CGs in our sample using the most precise 



velocities available from NED, the two galaxy median separators from Hickson et al. (1992), and 



total group H I masses from Green Bank Telescope, Arecibo Observatory, and Effelsberg 100 m 



Antenna single dish measurements by both Verdes-Montenegro et al. (2001) and Borthakur et al 



(2010 Subtypes based on the morphology of the H I gas, which we qualitatively assessed from 
VLA observations, are included as part of the proposed evolutionary sequence of compact groups 
( Konstantopoulos et al.|2010 ). Type A groups are those in which the neutral gas is confined to the 
individual group members, while Type B groups have H I gas distributed between the galaxies and 
not centered on any particular member (s). We first test the use of the H I to dynamical-mass ratio 
as a descriptor for group evolutionary state by comparing it to the group E/S0 galaxy fraction for 
the main galaxies. Using the H I mass to dynamical-mass ratio to characterize the evolutionary 
state is preferred over the E/S0 fraction because the mass ratio is a continuous distribution, while 
the E/S0 fraction values are discrete due to the small number of relatively massive members in each 



5 The H I masses of HCGs 16 and 62 are lower and upper limits, respectively. The lower limit of the mass in 
HCG 16 is due to the large angular size of the H I envelope, which extends beyond the Green Bank Telescope beam. 
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Fig. 5. — The H I gas mass ratio as a function of the main group member morphologies. The H I 
mass relative to the group dynamical mass decreases as the fraction of E/SO galaxies increases. The 
general trend of the data points indicates that CGs exhaust their H I reservoirs (i.e., Mm < 0.1% 
Mdyn) w h en approximately 50% of the main group members have E/SO morphologies. 



group. To quantify the relationship, we used the ASURV statistical package (Lavalley et al. 1992), 



which implements the methods presented in Isobe et al. (1986), to compute the Spearman rank 



correlation coefficient. This test measures how well the data are fit by a single monotonic function. 
From Figure [5j we find that the H I mass normalized by the group dynamical mass decreases with 
increasing E/SO galaxy fraction (67% probability from Spearman test). This result is expected if 
the H I to dynamical-mass ratio is indeed a tracer of the evolutionary state of the system. 

When we examine how Lx changes with H I mass ratio, we observe that the H I poor Type III 
HCGs 42, 59, 62, and 90 are X-ray brighter compared to more H I rich groups as expected if the 
gas has been processed by star formation or heated by the group potential. Of the Type II groups, 
we note that HCG 92 appears to be nearly as bright as the H I poor CGs, probably due to the 



shock front created by the high velocity intruder galaxy in the group (e.g., Trinchieri et al. 2003 
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Fig. 6. — Distribution of diffuse X-ray iuminosity in CGs as a function of the ratio of H I mass to 



dynamical-mass used in the evolutionary typing scheme defined by Johnson et al. (2007). Neutral 



gas masses are taken from Green Bank Telescope, Arecibo Observatory, and 100 m Effelsberg 



Antenna single dish measurements by Verdes-Montenegro et al. (2001) and Borthakur et al. (2010), 



while the H I morphologies are qualitatively assessed from VLA interferometric observations, which 
we lack for HCGs 59 and 62. Note that the values of the H I mass in HCGs 62 and 16 are upper and 
lower limits, respectively, and the plotted H I mass ratio corresponds to these limits. Futhermore, 



uncertainties in the H I masses from Verdes-Montenegro et al. (2001 ) are not available, therefore we 



only plot error bars in the H I mass ratio for data from Borthakur et al. (2010). After separating 



the groups by H I distribution subtype based on Konstantopoulos et al. (2010), there appears to 



be a distinction in the X-ray luminosity between the two populations. Due to the small number of 
groups in our sample, more data are required to concretely determine if this difference is real. 



O'Sullivan et al. 2009). The only H I rich Type I in our study is HCG 31, which contains diffuse 



X-ray emission linked to star formation activity. Figure [6] shows the X-ray luminosity as a function 
of the H I mass ratio. From the figure, we can see that compared to the hot gas in other Type III 
groups, HCGs 59 and 90 are 1-2 orders of magnitude fainter than the remaining H I -poor HCGs. 
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Fig. 7. — X-ray luminosity as a function of the specific star formation rate for the seven groups with 
SSFRs from Tzanavaris et al. (2010). The groups primarily fall into two classes: quiescent, X-ray 
brighter systems (solid box); and star- forming, X-ray fainter systems (dotted box). The boxes are 
only used to identify the groups in these two regimes, therefore the absolute positions and sizes 
of the boxes in the figure do not necessarily carry physical meaning. At the low and high SSFRs, 
the processes that give rise to the X-ray emission are dominated by the group potential and local 
means (e.g., superwinds), respectively. 



In the case of HCG 90, this could stem from the exclusion of galaxy A from the X-ray analysis due 
to its bright Seyfert 2 nucleus (see Appendix |A|), however it is unlikely that including the diffuse 
emission from 90A would increase the total group luminosity by an order of magnitude to bring it 
to the level of the brighter Type III groups. 

From Figure[6j we note that if we include the H I distribution subtypes of the groups, then there 
may be two distinct trends between the total diffuse X-ray luminosity and the H I evolutionary 
type. Specifically, the Type A groups (H I confined to galaxies) appear to be consistently less 
luminous in X-rays compared to Type B groups (H I in the IGM). Based on this tentative result, it 
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is reasonable to predict that HCG 62 should have any remaining H I dispersed throughout its IGM, 
while H I gas in HCG 59 should be mostly contained in and around its group members. However, 
we note that this result is preliminary and that more data are required to quantitatively assess the 
likelihood that the two H I gas morphologies are distinct populations in X-ray luminosity. If we 
assume these H I morphology classifications for HCGs 59 and 62 to increase the number of data 
points in each sample to the minimum required, a two sample K-S test gives a 95% probability that 
the H I subtypes are two distinct populations in X-ray luminosity. To examine how our assumption 
concerning the H I subtypes of HCGs 59 and 62 influenced the test, we perform a second K-S test 
in which we switch the H I morphology classifications for these groups (i.e., HCG 59 has intragroup 
H I, while 62 has H I only in the galaxies). This results in an 12% probability that the two H I 
morphologies represent distinguishable populations in Lx- Due to the small sample size, we cannot 
determine if this distinction is statistically significant, and further data are required to clarify this 
phenomenon. 

We note that the mix of galaxy morphologies in the groups complicates the comparison of the 



H I and X-ray gas. However, Verdes- Montenegro et al. (2001) and Borthakur et al. (2010) reported 



on the observed H I masses in these groups and the predicted H I masses based on the group member 



luminosities and morphologies using relations from Haynes & Giovanelli (1984) who examined H I 



in a sample of isolated galaxies. From this, we can compute H I deficiencies for the groups in our 
sample, where deficiency describes the ratio of the H I mass observed to that predicted (in contrast 
to our use of H I rich and poor, which describes the relative neutral gas mass normalized to the 



dynamical mass). Using the nomenclature of Borthakur et al. (2010) (i.e., "heavily" deficient, 



"slightly" deficient, and normal groups contain < 1/3, between 1/3 and 2/3, and > 2/3 of their 
predicted H I, respectively) we find that the X-ray fainter groups are either heavily deficient (7, 22, 
and 90) or normal (59) groups, while the brightest groups are either slightly deficient (16, 31, 62^j 
and (92) or normal (42) groups. 

Because the gas is not in hydrostatic equilibrium (as indicated by the diffuse X-ray morphology; 



see £4.1), and in many cases linked to the individual group members rather than the environment, 



we cannot use the temperature to calculate the group hot gas mass (such a procedure assumes 
hydrostatic equilibrium in three dimensions). However, if we did assume that the gas is in hy- 
drostatic equilibrium, then the narrow range of temperatures implies a similarly narrow range of 
masses ( Fabricant et al.||1980 ). The X-ray contours in Figure [T] make this unconvincing due to the 
extremely large number densities and/or metal abundances necessary to make the hot gas masses 
of groups such as HCG 16 similar to ones such as HCGs 90 or 62. If we instead hypothesize that the 
X-ray gas mass is proportional to the bolometric X-ray luminosity, i.e., that there exists a Lx~Mx 



relation for groups (as in Zhang et al. 2011 for clusters), then we can compare the relative X-ray 
luminosities of the groups as a proxy for their relative hot gas masses. That the X-ray fainter 



6 This group is likely heavily deficient as its H I mass is actually an upper limit and the ratio of the limit to the 
predicted value is 35% ( Verdes-Montenegro et al.|2001 1 . 
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groups for a given H I mass ratio are deficient in H I and that all of the CGs have similar X-ray 
temperatures suggests two scenarios: either (1) additional X-ray gas has too low surface brightness 
to detect in the available data, or (2) some fraction of it is missing from the groups. If the H I 
was not converted to hot gas, then where is the missing gas in compact groups? The least mas- 
sive groups in our sample have velocity dispersions that correspond to virial temperatures of only 
~0.08 keV, however if the gas is cooler than the virial theorem implies, then it may be in the form 
of a diffuse UV intragroup medium. Further data are required to fully explore this phenomenon. 



4.4. Diffuse X-ray Emission and Specific Star Formation Rates 

We also compare the X-ray emission of seven groups in our sample to the total group SSFRs 
for the main group members calculated from the UV+24 /jm fluxes measured by |Tzanavaris et al. 



(2010) with corrections from Tzanavaris (2012, private communication) in Figure [7J and find that 



there is a distinction in Lx for detected groups on either side of the SSFR gap (data are not 
available for HCGs 90 and 92). We exclude dwarf group members with measured SSFRs because 
HCG 31 is the only CG in our sample that has such measurements of its dwarf population, which 
increases the total group SSFR by several orders of magnitude due to the combination of relatively 
low to moderate SFRs and small stellar masses in star-forming dwarf galaxies. In our sample, 
CGs containing star-forming galaxies all share similar X-ray luminosities and are 1-2 orders of 
magnitude fainter than groups with low total SSFRs. This may indicate an "X-ray gap" analogous 
to the SSFR and mid-IR color gaps found by |Tzanavaris et al-pOlO] and |Walker et aTJ|2uT0"l |2"0~1~2 



We note that HCGs 7 and 22 have SSFRs that lie within the gap range of (0.3 - 1.8) x 10~ n yr" 1 , 
and neither has any detected diffuse X-ray emission. 

The brighter Lx values associated with the low SSFR groups are due to the advanced evolu- 
tionary stage of these environments. HCGs 42 and 62 have very low H I to dynamical-mass ratios, 
while simultaneously their relatively massive galaxies are entirely E/S0 types. From the evolution- 



ary scenario presented by Verdes- Montenegro et al. (2001 ), and the fact that compact environments 



are favorable to the tidal stripping of gas from galaxies (Freeland et al. 2009), we should expect 
that these older compact groups have removed the cool gas from the galaxies and heated it in the 
intragroup medium. Conversely, the star-forming groups HCGs 16, 31, and 59 have relatively faint 
X-ray emission associated with them. These groups also all exhibit only galaxy-linked emission, as 
expected from their relatively young evolutionary states. Finally, neither of the intermediate SSFR 
groups, HCGs 7 and 22, are detected by Chandra. From the figure, we find that there are two 
possible evolutionary scenarios with respect to the total group SSFR: (1) the groups move from 
the lower right portion of Figure [7] to the upper left, i.e., star-forming and X-ray fainter to quies- 
cent and X-ray brighter, though perhaps not monotonically; or (2) there exists at least one more 
evolutionary track in which star formation in CG galaxies declines while gas is not stripped from 
the disks and/or heated. In particular, studies of groups similar to HCGs 7 and 22 may provide 
further insight. 
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5. Summary 

We detect diffuse X-ray emission in seven of nine of the CGs in our sample with temperatures 
ranging from 0.6-0.72 keV and bolometric X-ray luminosities between 10 404 and 10 42 2 erg s _1 . 
The groups exhibit a wide range of velocity dispersions (56-343 km s _1 ), log 10 (Mn i)/ log 10 (Md yn ) 
(0.70-0.97), and morphological fractions from spiral-only groups to systems rich with E/S0 galaxies. 

Based on the hot gas morphologies, we find that the X-ray emission likely arises due to both 
local processes (i.e. star formation, nuclear activity, and tidal interactions) and global processes 
(i.e., heating by the group potential). In dynamically unevolved (i.e., low- mass, low velocity dis- 
persion) systems, the observable diffuse X-ray emission is dominated by local processes. The X-ray 
brighter groups (for a given H I mass ratio) have emission that stems from both an extended diffuse 
component (i.e., a true intragroup medium) and galaxy- linked emission, while emission detected in 
the X-ray fainter groups is only associated with the individual galaxies. 

HCGs 31, 42, 59, and 92 have X-ray luminosities in agreement with the predicted values from 
the Lx — T scaling relation from clusters, though the error in the temperatures of HCGs 31 and 59 
are large. Furthermore, HCGs 16 and 31 agree with the Lx—cr relation. In both Lx~T and Lx 
HCG 62 appears to lie within the scatter of the cluster data and is similar to fainter Lx clusters. 
The agreement between HCGs 16 and 31 with the Lx — o relation are likely coincidental because 
the scaling relation is predicted from the Virial theorem (i.e., Lx oc <7 4 ), and the X-ray emission 
from these two systems is clearly dominated by star formation rather than virilization. When the 
groups disagree with the cluster scaling relations, particularly in Lx — cr space, this indicates that 
the groups are not simply scaled-down analogs to galaxy clusters. Furthermore, given that the 
hot gas in the low-mass (i.e., low velocity dispersion) systems is found to be isolated to the group 
members rather than throughout the intragroup medium, we conclude that galaxy clusters are not 
a proper comparison class of objects for these groups (noting the possible exception of HCG 62). 

We also find that there may be a relationship between Lx and how the H I gas is distributed: 
preliminary evidence suggests that CGs with gas stripped from the galaxies are brighter in X-rays 
than groups with H I confined to the members, possibly due to strong multi-galaxy interactions 
that dispersed neutral gas into the intragroup medium and triggered star formation. However, 
the X-ray faintest groups are also more heavily deficient in H I implying that there may be some 
fraction of missing gas, possibly too cool to emit in X-rays or with too low surface brightness to 
detect. Finally, we note that groups dominated by local heating mechanisms have high UV+24 /im 
specific star formation rates, while groups with gas heated by the group potential have low SSFRs. 
The values of Lx between these two categories spans ~2 orders of magnitude and may indicate the 



presence of an 'X-ray gap" in CGs similar to the SSFR and mid-IR color gaps found by Tzanavaris 



et al. (2010) and Walker et al. (2010, 2012) 



The faintest Lx groups appear to be at very early stages in their evolution, perhaps coming 
together for the first time as is indicated by their low fractions of E/S0 galaxies. The influence 
of multi-galaxy interactions on liberating neutral gas from the galaxies and depositing it into the 
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intragroup medium early in the group lifetime (e.g., as seen in HCG 16) appears to have an effect 
on the ability of these groups to evolve into more cluster-like systems with respect to the hot gas 
distribution (e.g., HCG 62). 

5.1. Future Work 

Expanding the sample with appropriate observations of groups from the Chandra and XMM- 
Newton archives is the logical next step. This will give us a larger sample with which to study 
the relation between the hot gas and the evolution of group environment (e.g., how the X-ray 
luminosity varies with H I mass ratio). Furthermore, continuing to examine the differences between 
compact, loose, and fossil groups will demonstrate how efficient the compact environment processes 
gas. Inclusion of multiwavelength data will help to facilitate comparison of the hot gas to gas in 
cooler states (e.g., cold molecular). HCG 16 in particular is an interesting group worth more study; 
the presence of two potential M82-like superwinds in a single system presents an interesting case 
study for tidally induced star formation and how superwinds drive hot gas into the intragroup 
medium. 
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A. Notes on Individual Groups 

In addition to the unique aspects of the X-ray analysis and properties of the diffuse emission, 
we also list the center RA and Dec, shape, and dimensions of the extraction region for each group. 

HCG 7. — Diffuse emission in this group was not detected above the background in the Chandra 
data. The circular extraction region with radius 3'.9 was centered at RA = 00 h 39 m 23 s .9 and 
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Dec = +00° 52' 15 4". 

HCG 16. — The shallow observation of this group prevented detection of a true intragroup 



medium (i.e., the Lx = 2.4 x 10 40 erg s 1 IGM, corrected to our cosmology, found by Belsole et al. 



2003 with 45 ks of XMM data); however, hot gas associated with the individual galaxy members 
was detected. The extracted spectrum corresponds to the area surrounding galaxies A, B, C, and 
D, but does not include galaxy X, which is far removed from the group center and located far from 
the S3 aimpoint of the observation. The rectangular 7f2x3'7 extraction region was centered at 
RA = 02 h 09 m 33 s .O and Dec = -10° 09' 05 8" with PA = 36°. 

HCG 22. — Diffuse emission in this group was not detected above the background in the Chandra 
data. The rectangular 6?5x6'9 extraction region was centered at RA = 03 h 03 m 30 s .l and Dec = - 
15° 39' 27 3" with PA = 22°. We note that this region includes a background pair of galaxies to 
the southeast of the main group members, however any emission from this pair is negligible and 
does not adversely effect the upper limit on the X-ray luminosity of HCG 22. 

HCG 31. — The S/N of the detection is only marginally above the threshold required for a de- 
tection (S/N= 3), and therefore the properties of the hot gas in this system are poorly constrained. 
However, this is the first detection of diffuse X-ray emission in this group. The extraction region 
covers the massive group members and the southern tidal tail where the majority of star formation 



is occurring within the group (Gallagher et al. 2010). Based on the reservoir of H I gas in the 



group, Gallagher et al. (2010) predicted that conversion of 75% of the neutral gas into stellar mass 
over a 150 Myr episode of star formation would generate ~2 x 10 39 erg s" 1 of X -ray emission, 
however this is an order of magnitude below the observed value of Lx from the Chandra data. The 
elliptical 2^1x1^3 extraction region was centered at RA = 05 h 01 m 39 s .7 and Dec = -04° 16' 15 8" 
with PA = 44°. 

HCG lj.2. — We extracted spectra from both the region around galaxy A, where the diffuse 
X-ray emission is most readily apparent, and from a region containing all of the massive group 
members. However, the count rates obtained from both extraction regions were consistent within 
errors. Therefore, we find that the majority of the X-ray emission in the system is associated with 
galaxy A. To ensure that we do not exclude any extended emission, we compare our value with that 



derived from ROSAT PSPC data using the extended radius of the X-ray emission from Mulchaey 



& Zabludoff (1998). This radius corresponds to the distance at which the X-ray emission falls to 
20% of its peak value; in HCG 42, this radius is 8'. Comparison with the flux extracted from the 
ROSAT data shows that our value of the flux within our extraction region matches that obtained 
with an 8' radius. The circular extraction region with radius 3'A was centered at RA = 10 h 00 m 21 s .0 
and Dec = -19° 38' 48 5". 

HCG 59. — As for HCG 31, the X-ray emission from this group is extremely weak; however, the 
S/N is sufficient to classify this group as a detection. Due to the low S/N, the values derived from 
model fits to the extracted spectra are poorly constrained. Though it is unlikely that the value of 
the temperature could span two orders of magnitude (as indicated by the upper 90% confidence 
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error estimate), we include this error in subsequent figures for consistency. The 3^9x2' elliptical 
extraction region was centered at RA = ll h 48 m 26 s .5 and Dec = 12° 43' 10 2" with PA = 0.3°. 

HCG 62. — Previous work has found that there are cavities in the X-ray emission around 
HCG 62 (the result of AGN jets and lobes due to the Seyfert 2 nucleus in galaxy A) that lack high 



frequency radio emission, but do have powerful low frequency emission (Dong et al. 2010, Gitti 



et al.||2010[ |Giacintucci et al.||201l] |Q'Sullivan et al.||201l"| ). The extraction region does not include 
galaxy D, which is far to the south with respect to the other group members, and has negligible 
X-ray emission. During model fitting, we found that a model with a single MEKAL component was 
insufficient to properly fit the observed spectrum, therefore we model this group with two separate 



plasmas. We note that Mulchaey & Zabludoff (1998) found that the diffuse emission extends much 



farther from the group center than is evident in the Chandra data, and well beyond the FOV of 



the ACIS CCDs. Using the extended X-ray emission radius of 24'2 from Mulchaey & Zabludoff 



(1998), we found that our measured absorbed flux was a factor of 3.1 lower than that found from 



the ROSAT PSPC data and thus correct our measurements by this amount. Note that we apply 
this correction to both components of the emission in addition to the total luminosity, therefore 
the luminosities of the hot and cold component should be considered upper limits. In Figure [3j we 
separately plot both the cooler and hotter components. The circular extraction region with radius 
1^5 was centered at RA = 12 h 53 m 06 s .0 and Dec = -09° 12' 11 6". 

HCG 90. — The spectral extraction region for this group is centered on galaxies B, C, and D 
that are currently interacting with one another and are embedded within a halo of diffuse optical 
light (White et al. 2003). The region excludes the brightest member (90 A), which is located ~6'9 



(68 kpc) from the other three massive, interacting group members and contains a powerful Seyfert 2 
nucleus. This bright AGN complicates analysis of the diffuse emission with bright readout streaks 
and substantial pileup on the 13 CCD. The 5^3x4^1 elliptical extraction region was centered at 
RA = 22 h 02 m 04 s .5 and Dec = -31° 58' 51 9" with PA = 100°. 

HCG 92. — The known primary source of heating for the hot gas is a shock front caused by 
the high velocity intruder galaxy NGC 7318B as it moves at ~850 km s _1 through the intragroup 



medium (Pietsch et al. 1997). Numerous interactions have occurred in the group in the past 



500 Myr leading to tidal tails and debris (e.g., Fedotov et al. 2011 Hwang et al. 2012). These 



frequent interactions in the group likely caused gas to be stripped from the member galaxies and 



be deposited into the intragroup medium (Moles et al. 1997 Guillard et al. 2012). We note that 



the presence of the shock in the intragroup medium could have non-thermal X-ray emission that is 
not included in the MEKAL model fit. The circular extraction region with radius 2'.f> was centered 
at RA = 22 h 35 m 59 s .5 and Dec = 33° 58' 03 2". 
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